Dimensionality dependence of the wave function statistics at the Anderson transition 



O 

o 

(N 

> 
O 



CZ2 



X3 

o 



oo 

o 



O 



o 



X 



A. Mildenberger\ F. Evers^, and A. D. Mirlin^'^'* 
^Institut fur Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
^Institut fur Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 

(February 1, 2008) 

The statistics of critical wave functions at the Anderson transition in three and four dimensions 
are studied numerically. The distribution of the inverse participation ratios (IPR) Pq is shown to 
acquire a scale-invariant form in the limit of large system size. Multifractality spectra governing the 
scaling of the ensemble-averaged IPRs are determined. Conjectures concerning the IPR statistics 
and the multifractality at the Anderson transition in a high spatial dimensionality are formulated. 
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A disordered electronic system in d > 2 dimensions can 
be driven from the phase of extended states to that of lo- 
calized states by increasing the strength of disorder. This 
transition, bearing the name of Anderson, is character- 
ized by remarkably rich critical properties. In particular, 
the eigenfunctions at the critical point show strong fluc- 
tuations and represent multifractal distributions. 

These fluctuations can be quantitatively characterized 
by a set of inverse participation ratios (IPR), Pq = 
J d'^r |-!/'(r)p'. Using the renormalization group in d = 
2 + e dimensions (e ^ 1), Wegner found that the 
ensemble-averaged IPRs, (Pq), show at criticality an 
anomalous scaling with respect to the system size L, 
(Pq) cx where [| 



q{q-l)e + Oie'). 



(1) 



Equation (|l|) is written for the case of unbroken time- 
reversal symmetry corresponding to the orthogonal en- 
semble which we consider in the paper. According to 
Eq. (j^), the fractal dimensions Dq — fq/(q~l) are differ- 
ent from the spatial dimension d and depend on q, man- 
ifesting the multifractal character of the wave function 
intensity |'!/'(r)P [||. It is customary to characterize such 
a distribution by its singularity spectrum determined by 
the Legendre transform of f,, yielding 



/(«) 



id + e~af/Ae + 0{e^) 
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In the present paper we study how the wave func- 
tion statistics and, in particular, the multifractality spec- 
trum evolve with increasing spatial dimensionality, when 
the transition shifts from the weak-disorder range (as 
in d = 2 -f e) to the strong-disorder one. A numerical 
study of the Anderson transition in higher-dimensional 
systems is of special interest since no analytical results 
for the critical behavior in high d are available. In con- 
trast to conventional second-order phase transitions for 
which the mean-field treatment becomes valid above the 
upper critical dimension dc, so that the critical expo- 
nents are d-independent for d > dc, the usual mean-field 
approach fails in this case. The model has been solved on 
the Bethe lattice [|| which is believed to correspond to the 



limit d = oo. While for the critical exponent of the local- 
ization length the conventional mean- field value v = 1/2 
was obtained, a very unusual non-power-like critical be- 
havior of other quantities was found (namely, an expo- 
nential vanishing of the diffusion constant and a jump in 
the inverse participation ratios (Pq) with q = 2, 3, ... at 
the mobility edge). It was, however, shown in Q that 
these results are intimately related to the spatial struc- 
ture of the Bethe lattice. On this basis, it was argued 
that although the symmetry-breaking description of the 
transition in terms of an order-parameter function |^ is 
of general validity, the peculiar critical behavior found 
in is an artifact of the Bethe lattice and should take 
a power-law form for any finite d < oo. Assuming that 
the critical indices should match at d — > oo the Bethe 
lattice behavior (which corresponds to = for g > 2 
and to the value s = oo for the conductivity exponent), 
one then concludes that the upper critical dimension is 
dc = oo (see also ft). While available numerical re- 
sults for the value of and for the form of the critical 
level statistics Q in four dimensions (4D) are consistent 
with this conjecture, its rigorous justification and a sys- 
tematic analytical study of the transition for large d are 
still missing. 

To investigate the dimensionality dependence of the 
critical statistics of wave functions, we have performed 
numerical simulations of 3D and 4D tight-binding mod- 
els with periodic boundary conditions and a box distri- 
bution of site energies. The critical values of the disorder 
are known to be Wc ~ 16.5 for 3D @ and Wc ~ 35 
for 4D |^,||,|ll|. We have calculated wavefunctions with 
energy close to zero by diagonalizing the Hamiltonian 
using efficient numerical packages ||l^,|l3|. Thereby we 
could average over an ensemble that contained typically 
10'^ samples. From every sample 128 wavefunctions with 
energy close to zero have been taken into account. 

We begin by showing in Fig. 1 the evolution of the 
distribution 'P(lnP2) in 3D with the linear size L of the 
system from L = 8 to L = 80. It is clearly seen that 
at large L the distribution acquires a scale-independent 
limiting form and simply shifts along the x-axis with- 
out changing its shape. The scale-invariance of the IPR 
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distribution at criticality was conjectured in ||T^. This 
conjecture was questioned in [T5[ |, where numerical sim- 
ulations of the IPR distribution in 3D were performed, 
with the conclusion that the fractal dimension D2 is not 
a well defined quantity but rather shows strong fluctu- 
ations with root mean square deviation rms(I?2) 1- 
In terms of the IPR distribution the statement of 
would mean that the rms deviation (72 = [var(ln ^2)]^^ 
scales in the same way as (— lnP2) c>c InL in the limit 
of large L. However, a detailed analytical and numerical 
study of the wave function statistics in a family of crit- 
ical power-law random banded matrix (PRBM) ensem- 
bles has shown that the IPR distribution V{hiPq) is 
scale- invariant, corroborating the conjecture of It 
was argued that this is a general feature of the Anderson 
transition, and the conclusion of Rcf. was criticized 
as based on a not sufRciently careful numerical analysis 
of data obtained for too small systems. This expectation 
was supported by a very recent numerical study of the 
3D Anderson transition [pT| where a trend to saturation 
of the width of the distribution 'P(lnP2) with increasing 
length L was found. Our data (Fig. [|a) demonstrate the 
scale-invariance of the limiting distribution even more 
convincingly, since we have reached considerably larger 
values of L. The same conclusion can be drawn from the 
inset of Fig. |^ where the saturation of the rms deviation 
Cq{L) of the distribution ■p(lnPg) at L ^ cxd is demon- 
strated. Qualitatively similar results arc obtained for the 
IPR distribution in 4D (Fig. |l|b). Though a tendency to 
the saturation with increasing L is clear in this case as 
well, the full saturation has not been reached, in view of 
smaller linear sizes of the system as compared to 3D. 
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FIG. 1. IPR distribution in (a) 3D (system sizes L — 8, 
11, 16, 22, 32, 44, 64, 80) and (b) 4D (i = 8, 10, 12, 14, 16). 

To visualize the strength of the IPR fluctuations, we 
show in the main panel of Fig. ^ the values of the rms 
deviations Cg, which characterize the width of the distri- 
bution functions ■p(lnPg), extrapolated to L ^ 00. In 



d = 2 + e dimensions aq can be calculated analytically 
following iQ, with the result 



aq = Stt 026 q {q 



(3) 



where 02 — 0.00387 for the periodic boundary conditions. 
Here qc is the value of q corresponding to the root a_ 
of the singularity spectrum /(a), i.e. Qc = /'(a_) and 
f{a-) = 0. For qc the IPR distribution is dominated 
by its slowly decaying power-law "tail" 



P > ptyp 

^ ^q 



where Xq — qc/q This yields aq = q/qc 
dimensions, 

aq = (e/2)i/2q, q » (2/e)V2. 



(4) 

or in 2 -|- e 
(5) 

In addition to the numerical results for 3D and 4D, we 
present in Fig. |^ the result (||) of the e-expansion for 
d — 2.2. Furthermore, we show the small-g and large-g 
analytical asymptotics (||) and (||) with e = 1. Of course, 
the e-expansion is only justified parametrically for e <C 1. 
Nevertheless, we see that it still describes very reasonably 
the IPR fluctuations at the Anderson transition in 3D. 




FIG. 2. The rms deviation aq of InP^ extrapolated to 
L — > 00 in 3D (x) and 4D (stars). The dotted line is the an- 
alytical result (^) for e = 0.2; the full lines represent Eqs. (^) 
and (^) with e — 1. Inset: evolution of aq with L in 3D. The 
leading finite size correction of all data has the form with 
y = 0.25 ^ 0.5 for 3D and y = 0.1 0.4 in 4D. The numerical 
error in the extrapolated values of (7q(oo) is as big as 10% due 
to the uncertainty in y. 

Having demonstrated that the distributions of the 
IPRs Pq are scale-invariant at criticality, so that the frac- 
tal dimensions are well defined, we are prepared to an- 
alyze the form of the multifractal spectra. We use the 
numerical procedure described in Ref. [ p^ , where it was 
applied to the quantum Hall plateau transition. Specif- 
ically, we calculate the ensemble-averaged IPR (Pq) for 
different system sizes L, extract Tq and f{a) and extrap- 
olate to the thermodynamic limit L — > c» in order to 
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eliminate the finite-size corrections. We refer the reader 
to Ref. for more details of the procedure and for 
a discussion of its advantages as compared to the box- 
counting calculations performed in earlier publications. 




FIG. 3. Singularity spectrum /(a) in 3D (dashed) and 4D 
(full line). To illustrate the evolution of the spectrum from 
d = 2 to d = 4, analytical results for d = 2 + e are shown for 
e = 0.2 (dotted) and e = 0.01 (dot-dashed). Inset: compari- 
son between f{a) for 3D and Eq. (^) with e = 1 (solid). 

The obtained results for /(a) in 3D and 4D are shown 
in Fig. ^ To illustrate the evolution of the spectrum in 
the whole range from d — 2 to d = 4, the analytical re- 
sults (HI) for 2 -I- e dimensions with e = 0.2 and e — 0.01 
are also shown. In Fig. ^ the corresponding results for 
the fractal dimensions D{q) are presented. We see that 
with increasing dimensionality the singularity spectrum 
/(a) broadens. This is not surprising: with increasing d 
the transition moves further in the region of strong dis- 
order, implying stronger multifractality. What is much 
less obvious is that in the range a < 2 the /(a) curve 
shifts to the left with increasing d. This corresponds to 
the fact that the fractal exponents Dq with g > 1 de- 
crease with increasing d. In particular, for the exponent 
Z?2 determining the spatial dispersion of the diffusion 
coefficient at criticality we find D2 = 1.3 ± 0.05 in 3D 
and 1)2 = 0.9 ± 0.15 in 4D ||l^. It is worth mention- 
ing that the e-expansion (^ with e = 1 describes the 3D 
spectrum remarkably well (though with detectable devi- 
ations, see inset of Fig. ^). In particular, the position 
of the maximum, ao = 4.03 ± 0.05 is very close to its 
value ao = d + e implied by (||). As expected, in 4D 
the deviations from the parabolic shape are much more 
pronounced and ao = 6.5 ± 0.2 differs noticeably from 6. 

The obtained value of D2 = 1.3 ± 0.05 in 3D is consid- 
erably smaller than what was found in the earlier numer- 
ical studies 1 20 1 where the values in the range 1.4 1.8 
were reported. The reasons for this are as follows. In the 
earlier works the spectrum f{a) of individual eigenstates 
was studied via the box-counting procedure. While in 



the limit L ^ 00 the spectrum /(a) defined in this way 
should reproduce the part of the /(a) curve lying above 
the a;-axis, for a finite L the finite-size effects affect f{a) 
strongly, especially in the region close to the zero q;_ of 
f{a). As a result, such a method not only fails to yield 
the exponents Dq with q > but also leads to very 
large errors in determination of Dq with q smaller than 
but close to qc- In particular, in 3D and 4D we find that 
qc is close to 2, explaining large errors in earlier results 
for (7 = 2 in 3D. As to the 4D case, we are not aware 
of any previous studies of the wave function statistics at 
the Anderson transition. 




FIG. 4. Fractal dimensions Dq in 3D (dashed) and 4D (full 
line). Analytical results for d = 2 -\- e with e = 0.2 (dotted) 
and e = 0.01 (dot-dashed) are also shown. 

Let us discuss now what our findings imply for the 
high-d behavior of the critical wave function statistics. 
First of all, we notice that our numerical observation 
that Dq with q^l decrease with increasing dimensional- 
ity confirms the above-mentioned expectation that they 
should tend to zero in the limit d —>■ 00. In view of this, 
an analogy with the PRBM ensemble with the parameter 
6 ^ 1 is very instructive. This ensemble can be consid- 
ered as describing a ID chain with random long-range 
hopping whose rms amplitude is b/r where r is the dis- 
tance. The model is critical for arbitrary < 5 < 00 
and can be studied analytically in both limits 6^1 and 
6^1. The latter case is relevant to the issue under 
discussion, and we remind the reader of the key results 
p6[ on the wave function statistics at 6 <C 1. (The case 
6 3> 1 is to a large extent analogous to the Anderson 
transition in 2 -t- e dimensions with e <^ 1.) Specifically, 
the scale-invariant critical distribution V{Pq/Pq^'^) be- 
comes 6-independent in the limit 6^1. Furthermore, 
the exponents Tq with q > 1/2 are proportional to h in 
the small-6 limit. 



2bT„ 



f, = 2r(q~l/2)/7rl/2r(q-l). (6) 



Correspondingly, the singularity spectrum f{a) acquires 
for a 1 the form 
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/(«) = 2hF{a/2b) , 



(7) Hall-Systeme" der Deutschen Forschungsgemeinschaft. 



where F{A) is the Legendre transform of T{q), with the 
asymptotics F{A) ~ -1/nA as A ^ and F{A) ~ A/2 
as A — > DO. The smahness of the fractal exponents, 
fg ^ 1, reflects a very sparse structure of the eigen- 
states formed by resonances (spikes) with a hierarchy 
of distances between them, ri ^ r2 <g; . . ., such that 
lnr,+i/r, ^1/6 

After this reminder we return to the Anderson transi- 
tion in high dimensionality. The smallness of the fractal 
exponents, Dq ^ 1 at d ;5> 1 implies, in close similar- 
ity to the PRBM model, that the corresponding critical 
eigenstates have the resonance structure with a hierarchy 
of scales Inr^+i/ri ^ l/^d and with a dimensionality- 
dependent parameter bd satisfying &£( — * as d — > oo. 
(In contrast to the PRBM model, where the emergence 
of resonances was due to direct long-range hopping pro- 
cesses, now it should be determined by interference of all 
possible paths connecting two sites on a d-dimensional 
lattice.) The sparse resonance structure of eigenstates 
allows us to use the analogy with the PRBM model 
and to make the following conjectures: (i) The scale- 
invariant distribution V{Pql P^^) for q > 1/2 becomes 
essentially d-independent at d ^ 1. (ii) The fractal 
exponents take at d ^ 1 and for q > 1/2 the form 
fq{d) = 2bdTq, i.e. they depend on d only through the 
overall factor bd (satisfying 6^ ^ as d — > oo). Corre- 
spondingly, the singularity spectrum scales at large d as 
f{a) — 2bdF{a/2bd). Though we are fully aware of a 
speculative character of these conjectures at the present 
stage, we expect that they will stimulate further numer- 
ical and analytical work, which should finally resolve the 
challenge posed by the Anderson transition in high d. 

In conclusion, we have studied the statistics of criti- 
cal wave functions at the Anderson transition in 3D and 
4D. The distribution of the inverse participation ratio 
Pq was demonstrated to acquire a scale-invariant form in 
the limit of large system size. As a convenient measure of 
the strength of the IPR fluctuations, we have evaluated 
the rms deviation aq of In Pq and found a result match- 
ing well analytical predictions for 2 + e dimensions with 
e ^ 1. Calculating the ensemble-averaged IPR values, 
{Pq), we have determined the spectrum of fractal dimen- 
sions Tq = Dq{q — 1) and the singularity spectrum /(a) 
characterizing the multifractal properties of the wave 
functions. In particular, we have found D2 = 1.3 ± 0.05 
in 3D and D2 = 0.9±0.15 in 4D. More generally, our 
results indicate that the dimensions Dq with g ^ 1 de- 
crease with increasing spatial dimensionality d. On this 
basis, we have formulated two conjectures concerning the 
wave function statistics at criticality in the large-d limit. 
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